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ABSTRACT 

We present an analysis of multi-epoch low-resolution spectrophotometry, complemented by 
the light curves provided by massive photometric surveys spanning over 100 years, of the 
symbiotic binary LMC S63. We showed that it is an eclipsing binary with the orbital period 
of 1050“*. We also found evidence of outbursts in history of the white dwarf. If it was a Z-And 
type outburst, as is most likely, it would be a second such outburst recorded in the Magellanic 
Cloud symbiotic system. We confirmed that the red giant is enhanced in carbon, and estimated 
C/O ^ 1.2 by fitting a model atmosphere to the SALT spectrum. We also found bi-periodic 
pulsations of the red giant, and demonstrated that it is similar to other carbon variables with 
confirmed bi-periodicity. 
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1 INTRODUCTION 

Symbiotic stars are close binaries consisting of a late-type giant 
transferring material via wind or Roche-Lobe overflow to a much 
hotter compact companion surrounded by an ionized nebula. In 
case of a Mira as the red component a warm dust shell is present in 
the system (D-type symbiotic) in contrast to S-type with a normal 
giant (see e.g. Mikoiajewska 2012 for details). 

Currently there are 8 confirmed symbiotic stars in the LMC 
and 8 in the SMC (Belczyhski et al. 2000; Oliveira et al. 2013; 
Miszalski, Mikoiajewska & Udalski 2014). They are amongst the 
hottest and brightest discovered symbiotics (Mikoiajewska 2004). 
Additionally all of them contain asymptotic giant branch (AGB) 
stars, while both red giant branch (RGB) and AGB stars are found 
in the Galactic symbiotic systems (Mikoiajewska 2004). This could 
be caused by a selection effect and the situation could change as 
more results from ongoing surveys are published (e.g. Miszalski, 
Mikoiajewska & Udalski 2013; Miszalski & Mikoiajewska 2014; 
Rodriguez-Flores et al. 2014). 

Early observations of LMC S63 showed ~2 mag variability 
with a period of ~1020 days exhibiting irregularities around max- 
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imum (LMV 1485 light curve in Gaposhkin 1970)'. Based on this 
early light curve, Payne-Gaposchkin (1971) classified LMC S63 
as an R CrB variable. First spectra obtained by Feast & Webster 
(1974) showed a smooth continuum with high excitation emission 
lines, in particular H i, He i, [O iii] and broadened He ii, as expected 
for a symbiotic star. On further observations a late-type star con¬ 
tinuum was present confirming the symbiotic nature of the system 
(Allen 1980). Allen (1980) suggested that a carbon star could be the 
red giant component of LMC S63, which was later confirmed by 
Morgan (1992). Orbitally related variability in photometric obser¬ 
vations with a period of 1060“' was found by Mikoiajewska (2004). 

To study this star we carried out multi-epoch low-resolution 
spectrophotometry. The spectra were complemented by the light 
curves provided by massive photometric surveys. In Sect. 2 we de¬ 
scribe the data and their reductions. The photometric and spectro¬ 
scopic variability is analysed in Sect. 3. The physical parameters of 
the binary component are derived and discussed in Sect. 4. A brief 
summary of our results is given in Sect. 5. 


* Note that the period of 1200 days quoted by Feast & Webster (1974) is a 
printing error. 
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Table 1. The log of observations. The magnitudes are mean magnitudes 
predicted from the MACHO and OGLE V light curve (Fig. 3). 


Date 

MJD 

Phase 

Telescope 

Exposures [s] 

V 

14.11.2005 

53690 

0.62 

SAAO 1.9m 

2x1000, 1500 

15.8 

14.10.2006 

54022 

0.94 

SAAO 1.9m 

3x1000 

16.2 

21.10.2007 

54394 

0.29 

SAAO 1.9m 

3x1200 

15.8 

06.11.2008 

54776 

0.66 

SAAO 1.9m 

3x1200 

15.8 

18.11.2012 

56249 

0.06 

SALT 

2x1800 

16.0 

05.02.2013 

56328 

0.14 

VLT 

6x440 (UVB) 
6x350 (VIS) 
16x180 (NIR) 

15.9 


2 OBSERVATIONS 
2.1 Spectroscopy 

Our spectroscopic observations are summarized in Table 1. Low 
resolution spectra were made with Grating Spectrograph with the 
SITe (Scientific Imaging Technologies, INC.) CCD on the 1.9m 
telescope at the South African Astronomical Observatory (SAAO) 
during the period from 2005 to 2008. A long slit with a projected 
width of 1.5”, and grating #7 with 300 lines mm“* were used which 
resulted in a resolving power of RwlOOO. Every night at least one 
standard star was observed allowing for spectrophotometric cali¬ 
bration. Wavelength calibration was performed using CuAr refer¬ 
ence lamp spectra. The covered spectral range varied with every 
spectrum, but was approximately equal to 4000-7000 A. The stan¬ 
dard IRAF procedures were used to reduce the CCD observations 
and to perform flux calibration. Unfortunately, due to various rea¬ 
sons (like changes in sky conditions, slit losses, seeing) the flux 
calibration was affected by a gray shift. Therefore we scaled the 
spectra using known V magnitudes, and the cool giant features (see 
sec. 4.1) 

Another low resolution spectrum was obtained with the Robert 


Stobie Spectrograph (RSS; Burgh et al. 2003; Kobulnicky et al. 
2003) on the queue-scheduled Southern African Large Telescope 
(SALT; Buckley et al. 2006; O’Donoghue et al. 2006) under pro¬ 
gramme 2012-2-RSA_POL-001 (PI: Miszalski). Two 1800 s expo¬ 
sures were taken with the PG900 grating on 2012 November 18 to 
cover A ~4200-7275 A and A ~6155-9150 A. After basic calibra¬ 
tions were applied with the pysalt package (Crawford et al. 2010), 
we followed the same reduction process as outlined in Miszalski et 
al. (2014). The relative flux-calibrated spectrum was then normal¬ 
ized before scaling the spectrum to V = 16.0 mag, and the resultant 
spectrum is shown in Fig. 1. The emission line fluxes are listed in 
Table AL 

In addition, we used a publicly available spectrum obtained 
with X-shooter (Vernet et al. 2011) on the UT2 Very Large Tele¬ 
scope (VLT) at the European Southern Observatory under pro¬ 
gramme 090.D-0422(A). The spectrum was made simultaneously 
in three arms (UVB, VIS and NIR) which together covered the 
spectral range 3000-24000A. The slit widths for UVB, VIS and 
NIR arms were 1.0”, 0.9” and 0.9”, and resulting spectral resolu¬ 
tions R=4350, 7450 and 5300, respectively. 

We supplemented our own spectroscopy with 2 previously 
published spectra. The first one was obtained by D. Allen, in 1978, 
and analysed by Pereira (1995). The second one is from the spec¬ 
trophotometric atlas of symbiotic stars (Munari & Zwitter 2002, 
see therein for details). The emission line fluxes measured for these 
spectra are also included in Table AL In addition, we collected 
emission line fluxes published in Kafatos et al. (1983), Morgan 
(1992),Van Winckel, Duerbeck & Schwarz (1993) and Muerset, 
Schild & Vogel (1996). 

We measured emission lines in the UV on three Interna¬ 
tional Ultraviolet Explorer (lUE) spectra and one HST Faint Ob¬ 
ject Spectrograph (EOS) spectrum. Fluxes derived from three spec¬ 
tra were already presented in the literature (Vogel & Nussbaumer 
(1995) - EOS spectrum, Kafatos et al. (1983) - swpl6591, Pereira 
(1995) - swp25791), nevertheless since the flux calibration has been 
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improved over the years (see eg. Rodriguez-Pascual, Gonzalez- 
Riestra & Schartel 1999) we repeated the measurements. The re¬ 
sults are presented in Table A2. 

A typical error of the measured fluxes is 15% for the strongest 
and 30% for the weakest lines and lines in the UV range. This 
mostly represents uncertainty in selection of the level of the con¬ 
tinuum and flux calibration. The collected, previously published 
fluxes, had similar uncertainty. 

2.2 Photometry 

The photometric data in this study are from the OGLE (Udalski 
et al. 2008; Soszynski et al. 2011), MACHO (Alcock et al. 1991) 
and EROS-2 (Spano et al. 2011) projects. The typical error of mea¬ 
sured magnitudes in these surveys was 0.002 mag, 0.003 mag and 
0.05 mag respectively. The OGLE photometry is reduced to the 
Cousins V and I system. The MACHO data were transformed to 
the Cousins V and R using formulae from Lutz et al. (2010). The 
colour applied for transformation was a colour measured during 
individual nights. The EROS-2 data were taken simultaneously in 
the non-standard blue BE (420-720 nm) and red RE (620-920 nm) 
bands. The Harvard College Observatory archival plates covering 
the field including LMC S63, the same as those used by Gaposhkin 
(1970), are now digitalized within the framework of the DASCH 
(Digital Access to a Sky Century @ Harvard, Laycock et al. 2010) 
project. However, a comparison between the DASCH light curve 
and that published by Gaposhkin (1970) shows that some origi¬ 
nal plates must had been lost or perhaps not scanned. Therefore, 
we scanned Gaposhkin’s original light curve and scaled it to fit the 
DASCH B magnitudes from the same date. Error of magnitude esti¬ 
mates from photographic plates varied between 0.05 and 0.20 mag, 
which did not include uncertainty in the calibration process. 

Single photometric measurements for LMC S63 are also re¬ 
ported in the literature. These measurements are collected in Ta¬ 
ble A3. 

3 VARIABILITY 

3.1 Photometric variability 

The combined light curve of LMC S63 is presented in Fig. 2. This 
light curve is typical of eclipsing symbiotic stars which show deep 
and narrow eclipses during active phases, and more or less sinu¬ 
soidal light changes during quiescence (e.g. fig. 2 in Mikolajewska 
2003). 

The deep minima in the DASCH light curve can be reproduced 
with an ephemeris MJDmin = (27836.4±7.2) + (1049.7 ± 3.5) x£. 
Combining this ephemeris with the three minima visible in the MA¬ 
CHO and EROS light curves we found the orbital ephemeris to be 

MJDmin = (49885.6 ± 7.4) -t- (1050.0 ± 0.5) x E 

in agreement with the 1060'^ orbital period estimated by Mikola¬ 
jewska (2004) from MACHO light curves. Two of the eclipses 
in the period covered by the DASCH photometry are not de¬ 
tected; they were most probably missed due to their relatively short 
duration. The first Harvard observations were obtained in 1897 
(DASCH) and 1905 (Gaposhkin 1970), respectively, and in both 
LMC S63 has photographic magnitude /Kpg ~ 16. Whereas the 
DASCH point was taken near minimum {(j) = 0.07), Gaposhkin’s 
point is from the orbital phase of ~0.5, which strongly suggests that 
at that time the star was not yet in the active phase. 


In comparison to the most recent photometry (Fig. 3) the min¬ 
ima in the DASCH/Gaposhkin light curve are very narrow and have 
a much larger amplitude. This suggests that LMC S63 was in an ac¬ 
tive phase in 1930-54 when the DASCH/Gaposhkin light curve was 
obtained. 

Combining the OGLE, MACHO, and EROS data it is clear 
that the depth of minima during quiescence are larger in bluer fil¬ 
ters. This confirms that the variability is related to the hot compo¬ 
nent obscuration (Fig. 3). 

The OGLE light curves were also analysed by Angeloni et al. 
(2014), who found low amplitude changes with a 72-day period 
superposed on a general, slowly decreasing, trend. They attributed 
these periodic changes to red giant pulsations, and classified it as a 
LMC OSARG. However, they did not discuss the long-term trend 
because their baseline is too short. We can reveal with our longer 
time coverage that this decreasing trend is due to the orbital vari¬ 
ability. 

A power spectrum of the MACHO data after subtracting a si¬ 
nusoid with a period of 1050'* shows the most prominent peak at 
286'*. Unfortunately the MACHO data are not adequate for search¬ 
ing for 72“* variations because of the low time-resolution and the 
amplitude of the variability is similar to the accuracy of the mea¬ 
surements. EROS observations cover a period of time which par¬ 
tially overlaps the time period covered by MACHO. The EROS ob¬ 
servations confirm the 286“* pulsations (Fig. 3). Episodically the 72“* 
oscillation seems to also be detectable. Presumably during this time 
observations were carried out during the best observing conditions 
which allowed for detection of variability with such small ampli¬ 
tude. In the OGLE data the 286“* modulation is not found (Fig. 3). 
The cause of this could be the fact that the OGLE data mainly con¬ 
sists of a separate ~300'* datasets so it is not sufficient for searching 
for changes with such long periods. We conclude that in the time 
period covered by MACHO and EROS data (1993-2003) both the 
72'* and the 286“* variability was present but the 286“* variability 
most probably stopped before the OGLE observations (2010-2013). 
This would be in agreement with the fact that > 200“* variations su¬ 
perimposed to the variability on a time scale of tens of days could 
be episodic in a small-amplitude red variables, as was observed eg. 
in W Boo (Percy & Desjardins 1996). 

The average K = 11.3 locates LMC S63 above the First Giant 
Branch (FGB) tip, whereas the 286'* and 72“* periodicities are con¬ 
sistent with the Mira sequence and the B sequences, respectively, 
in the K log P plane (Wood 2000). This suggests that the longer pe¬ 
riod may represent the fundamental mode and the shorter one the 
first or a higher overtone. Based on studies of nearly 300 Galactic 
carbon giants, Bergeat, Knapik &. Rutily (2002a) discussed their 
physical parameters. Using the average < J — K >o=L14 and 
< K >0=11.3 for LMC S63 (Table Al), and adapting their calibra¬ 
tion of the bolometric correction BC(K) vs J-K colour we estimate 
Mboi=-4.4. Adopting the effective temperature Tetr = 3500 K con¬ 
sistent with the best fit to the optical spectrum and the IR SED (see 
section 4.1) we derive the radius Rg=184RQ, which is in very good 
agreement with the value resulting from our model fit to the optical 
spectrum, although both values could be affected by a systematic 
error. The red giant of LMC S63 then lies on the evolutionary track 
of a low metallicity (Z=0.008) 1 .9 Mq star in the HR diagram given 
in Fig. 9 of Bergeat et al. (2002a). 

The pulsation properties, and in particular the bi-periodicity of 
the Galactic carbon-rich variables was studied by Bergeat, Knapik 
& Rutily (2002b), who gave relations between their periods and 
physical parameters. In particular, the two periodicities found in 
S63 locate it in the period-radius diagram (Fig. 2 of Bergeat 
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Year 

1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010 



Figure 2. The DASCH and Gaposhkin (1970) (black and gray circles, respectively) light curve of LMC S63. The variability is caused by the eclipses (see 
text). The vertical lines indicate eclipse times according to our ephemeris. The an'ows indicate the dates of spech'oscopic observations discussed in the paper. 


et al. 2002b) among other carbon variables with confirmed bi¬ 
periodicity. However, its radius and the 286“* period is consistent 
with the fundamental mode only if the pulsation mass is ~ 0.8 Mq, 
much lower than the initial mass given by the position of S63 in 
the HR diagram. The ratio of the two periods, PolPi ~ 4 is much 
higher than the average value of ~ 2.24 found by Bergeat et al. 
(2002b), and the mean theoretical value, > 2 for fundamental to 
first overtone modes. 


3.2 Spectroscopic behaviour 

We present measured fluxes of the emission lines in our spectra in 
Table Al. Phase plots of the fluxes of selected emission lines are 
presented in Fig. 4. The fluxes of the H i and He i lines are lowest 
during the photometric minimum. Such behaviour is expected for 
orbitally related modulation and is caused by partial obscuration 
of the line forming region. Similar variability of the emission lines 
was observed, e.g., in Cl Cyg (Kenyon et al. 1991), and many other 
S-type symbiotic stars. Such orbitally related changes should be 
also observed in the Balmer UV continuum. In fact, the Galaxy 
Evolution Explorer (GALEX) (Martin et al. 2005; Morrissey et al. 
2007) observations of LMC S63 resulted in F(NUV)=647.59-b/- 
13.48 mJy on MJD 55803 ((/> = 0.700), and F(NUV)=553.52 -h/- 
7.92 mJy on MJD 55903 {(p = 0.796), and they seem to reflect 
mostly the expected changes in the Balmer continuum. 

The He 1 5876/6678 and 7065/6678 line ratios increase near 
the orbital phase (p = 0.0 (hot component obscuration; see Fig. 4). 
The reason for this is unclear but it is worth mentioning that these 
two lines are more likely enhanced due to a collisional excitation 
than the He 1 6678 line (see eg. Clegg 1987). The change in ratios 
of the lines are similar to those expected for changes associated 
with modulations of electron density in symbiotic binaries (Proga, 
Mikolajewska & Kenyon 1994). This is consistent with the Allen 
(1980) discussion about the fact that the high ratio of the Hei lines 


to the Balmer lines derived from his spectrum {(p = 0.95) indicates 
unusually high He abundance in the nebular region or collisional 
excitation of He i. 

He II 4686 shows no orbitally related modulations in contrast 
to the other lines (Fig. 5), except for the spectrum taken during the 
1978 eclipse (cp = 0.97) when it was only marginally detectable. It 
is possible that the eclipses in the He ii lines are very narrow or the 
line shows eclipse behaviour depending on its activity stage. Sim¬ 
ilar behaviour is observed in many other symbiotic stars, which 
show deep and narrow eclipses in Hen4686 only during some 
phases of their outbursts and nearly constant flux in quiescence 
(e.g. Cl Cyg; Kenyon et al. 1991). Moreover the line shows slow 
rise over the observed period (Fig. 6). Such behaviour can be ac¬ 
counted for by changing conditions and stratification effects in line 
forming region(s). 

The only forbidden lines present in our spectra are those of 
[Oiii]. Whereas [Oiii] 5007 is visible during the whole observed 
period [Oiii] 4363 was detected only in 1978, 1997 and 2006. Both 
lines show a gradual decline (Fig. 6). At the same time Hen4686 
shows a gradual increase. Similar behaviour was observed in Cl 
Cyg by Kenyon et al. (1991) during late decline from its major 
1975 outburst. 

We recall here that the first spectra of LMC S63 taken by 
Feast & Webster (1974) in 1972/73 (i^ ~ 0.2) revealed a broadened 
He II 4686 line with FWHM=:300km/s, which suggests formation 
near to the hot component, presumably in its wind. The line was 
even broader, with FWHM of 770 km/s, in late 1993 (Muerset et al. 
1996). The broadening of the He ii 4686 with respect to other nebu¬ 
lar lines, including H i, and [O iii] was also noticed by Allen (1980) 
and Kafatos et al. (1983). On the VLT spectrum the Hen4686 line 
is not so broad, its FWHM is only 134 km/s. Broadening of emis¬ 
sion lines is a common characteristics of symbiotic outbursts, and 
it is usually attributed to increased wind and expansion of the hot 
component atmosphere. LMC S63 very likely underwent an opti- 
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Figure 3. Left: Phase plot of the MACHO V & R data, OGLE V & I data (black and gray circles respectively), and EROS B & R data with the orbital 
ephemeris. Bottom panel: phase plot of the m^j photometric data from Gaposhkin (1970) and DASCH (gray and black circles, respectively) before and after 
MJD=25000 (open and filled circles respectively). Right: Phase plot with P=286'^ of MACHO V & R data and OGLE I data after subtraction of orbitally 
related variability. 


cal outburst in the 1930/50s. However, it is hard to say whether the 
presently observed secular changes are related to that event or an¬ 
other outburst had happened more recently. Nevertheless there was 
no additional activity of the system recorded in the archive data, 
although there is a gap from 1954-75 (see Table A3). 


4 PARAMETERS OE THE COMPONENTS 
4.1 The cool giant 

The cool component is a carbon star, which was first suggested 
by Allen (1980) and confirmed by Morgan (1992). Muerset et 
al. (1996) estimated its spectral type as type C2.1J using the 
monochromatic flux ratios defined by Cohen (1979) based on the 
classification scheme proposed by Yamashita (1972). However, the 
only attempts to determine the chemical abundances of this sys¬ 
tem were carried out for the nebula (Vogel & Nussbaumer 1995, 
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Figure 4. Phase plot of fluxes of selected emission lines. The fluxes are 
highest ai'ound phase 0.5 as expected for orbital modulation. However, the 
He I line ratios show a rise near phase 0 (the hot component behind the red 
giant). References: ♦ - this study, □ - Kafatos et al. (1983), ■ - Morgan 
(1992), X - Van Winckel et al. (1993), + - Muerset et al. (1996). 
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Figure 5. Comparison of orbitally related variability of He ii 4686 and H^ 
in the spectra taken in 2005-2012. 
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Figure 6. The gradual decrease of the [O iii] and C iii] line strength and in¬ 
crease of He II 4686 line. The O iii] line shows gradual decline superimposed 
to the orbital variability. References: • - this study, H - Kafatos et al. (1983), 
■ - Morgan (1992), -I- - Muerset et al. (1996). 
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Pereira 1995) and resulted in a carbon abundance too low for a car¬ 
bon enriched star. To determine the carbon abundance in the gi¬ 
ant we attempted fitting synthetic spectra to the observed spectrum. 
Synthetic spectra were calculated using Castelli & Kurucz (2003) 
models and the SPECTRUM code (Gray & Corbally 1994) with 
help of an updated Plez (1998) TiO lines list from the authors web¬ 
site^. The adapted solar abundances are from Grevesse & Sauval 
(1998). 

The input effective temperature, Tes, was estimated from the 
relation T^s = 7070I[(J - K) + 0.88] derived by Bessell, Wood & 
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Figure 7. Scaled model atmosphere spectra (black lines) with increasing 
carbon abundance compared to the SALT spectrum (red lines). 


Evans (1983) for K, M and C giants, for which J - K - 1.17 (Ta¬ 
ble Al) resulted in = 3450K. To compute the synthetic spec¬ 
trum we also assumed the mean LMC metallicity [M/H] = -0.5, 
and logg = 0 (appropriate for a low mass star with a radius of 
~ 200Rq, e.g. Muerset et al. 1996). 

After calculating a grid of models we performed a least 
squares fit of a scaled model to the SALT spectrum in the spec¬ 
tral range 5000A-6500A. As a result we obtained a best ht for 
[C/M] = 0.35 ± 0.02 dex and T,g = 3500 ± 250K (see Fig. 7). 
The scaling factor, {R^jcty^ = 7.48 x 10“^', corresponds to the red 
giant radius R^ = 190 Rq at the distance to the LMC, d = 49.97 
kpc (Pietrzyhski et al. 2013). 

The fit should not be affected by the nebular emission because 
the SALT spectrum was taken near to the minimum i(j> - 0.06). 
Moreover, the synthetic photometry of the model fits well to both 
optical and IR magnitudes confirming that the nebular contribu¬ 
tion is indeed negligible. When we assume the same oxygen abun¬ 
dance as the metallicity, and adopting the solar C/O=0.55 (Asplund 
et al. 2009), we receive C/O - 1.2, a clearly higher value than 
the one calculated for the nebula by Vogel & Nussbaumer (1995); 
C/O = 0.27 and C/O = 0.68 from 1982 and 1994 spectrum, re¬ 
spectively). Our result supports their suggestion that the C/O ra¬ 
tio in the nebula represents mostly the abundance of the hot com¬ 
ponent and possibly consists of material ejected during the active 
phase, which is consistent with the forbidden line variability. The 
increase of the C/O ratio in the nebula between 1982 and 1994 can 
be explained by the fact that the material in the nebula is slowly 
becoming dominated by material from the red giant wind. 

LMC S63 was included in the recent study of SAGE AGB can¬ 
didates in the LMC by Riebel et al. (2012), who determined T^ff 
of the cold component and the dust mass-loss rate by fitting model 


Table 2. Parameters of the cold component from fitting GRAMS C-rich 
models to IR data. The values from Bayesian analysis are presented. 



Best value 

Probability 

Teff 

3700 

54.60 % 

L 

9120 Lo 

54.60% 

logg 

0 

100.00 % 

Mass 

2 Mo 

100.00 % 

C/O 

2 

52.51 % 


2.53 X 10“" Mo/yr 

35.75 % 

Ein 

3 R-star 

51.78 % 

rii.s 

0.006 

32.30 % 


to spectral energy distribution based on 12 photometric bands from 
optical to mid-infared wavelengths. Unfortunately, they misclassi- 
fied LMC S63 as O-rich star, which could affect their fitting result. 
We repeated their procedure by fitting the GRAMS C-rich grid of 
models of dust shells around red supergiant and AGB stars (Srini- 
vasan, Sargent & Meixner 2011) to the same mid-infared photom¬ 
etry as Riebel et al. (2012), supplemented with data from WISE 
(Cutri et al. 2012) and DENIS (Cioni et al. 2000) catalogs. Using 
VOSA (Bayo et al. 2008) we performed a Bayesian analysis of the 
data. We assumed the distance of LMC d = 49.97 kpc (Pietrzynski 
et al. 2013) and E(B-V)=0.05 (Muerset et al. 1991). As a result, 
apart from parameters of the red giant, we estimated such parame¬ 
ters as the dust mass-loss rate Mdust, inner radius of the dust shell 
Rin and optical depth at 11.3 microns rii ,3 (see Table 2). Note that 
the grid of models available to us had a solar metallicity and high¬ 
est possible log g = 0, hence this results should be accepted with 
caution. Particularly the photospheric ratio of C/O-2 is higher than 
the value resulting from the optical spectra fitting, and it could be 
affected by the fact that LMC S63 has lower oxygen abundance 
than the one used in models. 


4.2 The hot component 

The published estimates of the hot component parameters were 
based on fitting a blackbody to the UV spectra in addition to 
the UV and optical Hen, and H/f emission line fluxes. Given the 
eclipsing nature of LMC S63, the only reliable estimates are those 
based on data collected out of the eclipse, and preferably near to 
phase ~ 0.5 when the hot component passes in the front of the 
red giant. In particular, Vogel & Nussbaumer (1995) derived the 
hot component temperature, T^, and luminosity, Li,, by fitting the 
He II 1640 emission line flux and the adjacent continuum. Their val¬ 
ues, 7’h=80-115kK and 80-90 kK, and L,, = 8200-11 500 L© and 
5000-6200 L© obtained from the 1982 lUE and 1994 lUE/HST 
data, respectively, are in fairly good agreement with other estimates 
(Kafatos et al. 1983; Muerset et al. 1991; Meier et al. 1994) based 
on the same lUE data. 

The maximum brightness of LMC S63, observed during the 
1930-1950 optical outburst, and corrected for a reddening E{B — 
V) = 0.05 (e.g. Muerset et al. 1991) consistent with the foreground 
reddening towards the LMC, is nipg 13.2. This corresponds to 
an absolute magnitude Mpg -5.3, which assuming that most of 
the hot component continuum emission is shifted to the optical (a 
lower limit if not), results in Li,,ma.x = 10 900Lo. 

In the 2000’s the hot component reached quiescence (as in¬ 
dicated by the optical brigthness, and the changes in the emission 
lines), and we can make use of the emission line fluxes listed in Ta¬ 
ble AL In particular, the lower limit on the hot component temper- 
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ature can be set via the relation between the highest observed ion¬ 
ization potential (IP) and the hot component temperature, namely 
r[K] 10^ IP[keV], proposed by Muerset & Nussbaumer (1994). 
On both our SALT and VLT spectra the highest detected IP was 
99 keV for Fe^®. An upper limit for T ~ 120kK was estimated from 
He It 4686, He 1 5876 and H^S emission line ratios observed outside 
of the eclipses (2005, 2007 and 2008 data in Table Al) assuming 
case B recombination and a black-body continuum (lijima 1981). 
Equations 6 and 7 of Mikolajewska, Acker & Stenholm (1997) then 
give L(Heii4686 - 1300 Lq, and L(H/}) == 1000 Lq. The temper¬ 
ature and luminosity of the hot component of LMC S63 is then 
similar to those of other queiscent symbiotic stars (e.g. Muerset et 
al. 1991, Mikolajewska et al. 1997). 

The high temperature, T > 100 kK, and luminosity, L > 
1000 Lq, indicates that LMC S63 could be a detectable soft X- 
ray source. Indeed, LMC S63 was detected at the 5cr level by the 
ROSAT PSPC survey in September 7-December 17, 1990 (0 = 
0.35 - 0.45; Bickert, Stencel & Luthardt 1993). The flux corrected 
for vignetting was 0.0112 counts/sec which is ~140 times less than 
the flux of AG Dra observed at quiescence on Nov 27-Dec 7 (Bick¬ 
ert et al. 1993), at orbital phase <p ~ 0.7 and L;, ~ 1000 (d/2.5 
kpc)^ (Mikolajewska et al. 1995). After scaling the AG Dra flux to 
the distance of the LMC (Pietrzyhski et al. 2013) its flux should be 
equal to ~0.00385 counts/sec, ~3 times fainter than LMC S63. We 
recall here that the LMC S63 hot component luminosity was ~ 5x 
higher in 1994 than that for AG Dra. 

The presence of broad He it 4686 in the spectra taken in 1974- 
1993, as well as the presence of broad high excitation lines in the 
HST spectrum taken in 1994 can be accounted for by the pres¬ 
ence of a fast wind from the hot component (Vogel & Nussbaumer 
1995). The weakness of the He it line during the 1978 eclipse of the 
hot component compared to that on the spectrum taken 2 years later 
by Kafatos et al. (1983) (see Fig. 6) additionally supports this inter¬ 
pretation. However, the He ii is relatively narrow on the VLT spec¬ 
trum taken in 2013, which indicates that the wind has ceased by this 
time, and the line does not show eclipses anymore. The change in 
the He ii characteristics coincides with a decrease of the hot compo¬ 
nent luminosity by almost one order of magnitude, which suggests 
that LMC S63 has recovered to a quiescent state. Unfortunately the 
~ 20-year gap in photometric observations does not allow us to as¬ 
certain whether the spectroscopic evolution was due to continued 
decline from the 1930-54 optical outburst. Nevertheless, we note 
that the hot component luminosity in 1982 and 1994 was compara¬ 
ble to that derived for the optical maximum. 

It is also not easy to deduce the nature of the outburst, whether 
this was a symbiotic nova or a Z And-type outburst. The light curve 
(Fig. 2) shows at least two maxima which is more typical for Z 
And-type outbursts (see, e.g., the light curves of Z And, AR Pav 
and other Z And-type systems in the AAVSO database). Some of 
the Z And-type systems also show very prolonged active phase 
(e.g. AR Pav has been active for more than a decade) or a series 
of outbursts (e.g. Z And). The outburst amplitude, Ampg ~ 2.5 
is also rather low compared to other symbiotic novae with their 
typical optical amplitudes exceeding ~ 5 mag (e.g. PU Vul - see 
fig. 8 of Kato, Mikolajewska & Hachisu 2012). The pre-outburst 
photographic magnitude, < Mpg 16 ± 0.3, is close to the B 
mag observed in the 1980-90s, which suggests that the object was 
brighter than expected for a C-rich red giant, and that the hot com¬ 
ponent and symbiotic nebula considerably contributed to the ob¬ 
served brightness. All these speak in favour of Z And-type outburst, 
which would make LMC S63 the second Z-And-type Magellanic 
symbiotic star after LIN 9 (Miszalski & Mikolajewska 2014), al¬ 


though we cannot completely exclude a weak nova outburst on a 
hot white dwarf. 


5 CONCLUSIONS 

We discussed for the first time orbital variability and outburst be¬ 
haviour based on photometry spanning more then 100 years, for 
Magellanic symbiotic system LMC S63. We showed that it is an 
eclipsing system with a orbital period of Porb = 1050.5 ± 0.5“*. 
We also found orbitally related changes in the emission line fluxes. 
The DASCH photographic light curve revealed a ~ 2.5-magnitude 
outburst which started in 1930. The star had remained active in 
the optical for at least 14 years. A broadened He ii 4686 with re¬ 
spect to other nebular lines revealed by the first spectra taken in the 
1970s, as well as all other spectra taken in 1980s and 1990s, sug¬ 
gests that either LMC S63 had undergone another outburst or the 
first outburst continued, and the system reached quiescence only in 
the 2000s. 

We demonstrated that the red giant is enhanced in carbon, and 
derived C/0 - 1.2 by fitting a model atmosphere to the SALT spec¬ 
trum taken near the optical minimum. We also discussed the pul¬ 
sation characteristics of the red giant, and compared them to other 
carbon variables. In particular, the two periodicities found locate 
the giant component of LMC S63 in the period-radius diagram 
among other carbon giants with confirmed bi-periodicity. 

Some of these conclusions results could not have been made 
by previous studies based on data spanning only a few years (e.g. 
Mikolajewska 2004, Angeloni et al. 2014), which also emphasises 
the need to look at data over a very long baseline to fully under¬ 
stand these systems. Because of its complex variability, LMC S63 
deserves to be monitored in the future. In particular, another out¬ 
burst would be among the most interesting thing to observe. 
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APPENDIX A: 


LMC S63: a deeply eclipsing symbiotic star 


Table Al. The observed wavelength of emission lines from the VLT spectrum and emission line fluxes in all of our spectra. 



Date 

15.02.1978 

17.10.1994 

14.11.2005 

14.10.2006 

21.10.2007 

06.11.2008 

18.11.2012 

05.02.2013 


MJD 

43555 

49643 

53690 

54022 

54394 

54776 

56249 

56328 


Orbital phase 

0.97 

0.77 

0.62 

0.94 

0.29 

0.66 

0.06 

0.14 

‘^ohs [A] 

ID 




Flux [10^'^' 

erg s“^cm“^] 




3124.6 

Oiii 3121.7 








0.23 

3135.6 

O III 3132.9 








0.95 

3205.6 

Hen 3203.1 








0.56 

3315.5 

Oiii 3312.3 








0.17 

3343.8 

Oiii 3340.7 








0.25 

3348.7 

[Nev] 3345.8 








0.37 

3428.7 

[Nev] 3425.9 








1.00 

3431.0 

Oiii 3428.7 








0.36 

3447.2 

Oiii 3444.1 








0.42 

3589.6 

[Fe vii] 3587.2 








0.08 

3706.6 

H 16 








0.02 

3724.8 

H 14 








0.03 

3737.2 

H 13 

0.10 







0.04 

3753.0 

H 12 

0.09 







0.05 

3758.2 

[Fev] 3755.7 

0.11 







0.08 

3762.4 

[Fevii] 3759.9 








0.12 

3763.4 

Oiii 3759.9 








0.21 

3773.6 

H 11 








0.08 

3800.8 

H 10 

0.15 







0.11 

3838.4 

H 9 

0.27 







0.25 

3872.2 

[Neiii] 3868.8 

0.42 

0.87 






0.07 

3892.2 

Hg.Hei 

0.50 

1.7 

0.53 





0.30 

3930.4 

He I 3926.5 








0.04 

3973.2 

H, 

0.46 

1.2 


0.64 




0.44 

4030.0 

He 1 4026.2 

0.23 







0.08 

4074.2 

[Fev] 4071.2 

0.13 







0.04 

4101.4 

N III 4097.3 

0.16 







0.14 

4105.0 

Hi 

1.1 

2.1 

4.0 

1.1 


1.70 


0.69 

4147.8 

Hei 4143.8 








0.08 

4203.8 

N III 4200.0 








0.06 

4343.8 

Hy 

1.7 

3.0 

6.1 

1.5 

1.6 

2.7 

1.6 

1.2 

4367.1 

[Oiii] 4363.2 

0.60 

0.54 


0.25 




0.06 

4392.2 

He 1 4387.9 

0.11 






0.13 

0.16 

4476.0 

He 1 4471.5 

0.64 



0.32 



0.29 

0.15 

4545.8 

He II 4541.6 







0.21 

0.12 

4638.6 

N III 4634.2 








0.08 

4643.4 

[Fe II] 4639.7 








0.05 

4645.2 

Niii 4640.6 

0.35“ 






0.18“ 

0.15 

4646.4 

N III 4641.9 








0.07 

4652.0 

C III 4647.4 

0.28'’ 






0.08'’ 

0.12 

4654.6 

Cm 4650.2 








0.09 

4689.8 

Hen 4685.7 

0.53 

2.7 

3.7 

3.1 

2.2 

2.6 

3.7 

2.9 

4718.0 

He 1 4613.2 

0.15 







0.09 

4865.0 

Up 

4.6 

9.3 

18.0 

5.0 

6.9 

8.3 

4.0 

2.9 

4926.8 

He 1 4921.9 

0.23 






0.26 

0.24 

5012.9 

[Oiii] 5006.9 

1.5 

1.3 

0.33 

0.60 


0.16 

0.15 

0.09 

5020.8 

He 1 5015.7 

0.37 






0.26 

0.23 

5416.7 

Hen 5411.5 







0.24 

0.24 

5881.6 

He 1 5875.6 

2.7 

2.0 

2.8 

3.2 

2.9 

1.8 

1.8 

1.1 

6091.8 

[Fevii] 6085.5 








0.16 

6094.4 

[Cav] 6086.4 








0.15 

6567.4 


29.6 

46.5 

72.0 

35.4 

66.0 

52.0 

22.3 

23.7 

6684.5 

He 1 6678.2 

1.1 

1.8 

2.1 

1.4 

2.9 

2.0 

0.85 

0.50 

7072.1 

He 1 7065.2 

2.6 

2.1 

2.6 

2.2 

3.0 

2.0 

2.0 

1.4 

10132.5 

Hen 10123.6 








0.68 

10841.3 

He I 10830.3 








3.3 


Notes: ^ blend of Niiid 4640.6, and /14641.9; ^ blend of Cm 4 4647.4, and /14650.2. 
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Table A2. UV emission line fluxes in units of 10 erg cm ^ s * 


Spectrum 

Date 

MJD 

Orbital phase 

swp16591 
21.03.1982 

45049 

0.38 

swp25791 

27.04.1985 

46182 

0.47 

FOS 

19.09.1994 

49614 

0.76 

swp53165 

23.12.1994 

49709 

0.85 

Nv 1240 

1.6 

1.3 

2.8 


Niv] 1486 

1.5 


1.4 


Civ 1554 

6.0 


7.9 

5.0 

Hen 1640 

2.4 

2.9 

2.4 

1.6 

Oiii] 1663 

1.5 

2.8 

0.8 


Niv 1720 

0.6 


0.6 


Niii] 1750 

1.1 


0.6 


Sim] 1893 

0.3 


0.1 


Cm] 1909 

4.2 

2.2 

1.6 

1.0 


This paper has been typeset from a TgX/ ETgX file prepared by the 
author. 
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Table A3. Photometry of LMC S63 from the literature. 


MJD 

Orbital 

phase 

B 

V 

R 

I 

J 

H 

K 

Reference 

42421 

0.89 

16.08-' 


14.62 

14.05 




1 

44284 

0.67 





12.48 

11.60 

11.33 

2 

44565'" 

0.93 

16.3 


14.1 





3 

44934'" 

0.27 

15.83 


14.29 

13.61 




4 

45049 

0.39 


15.4 






2 

45153'" 

0.48 

15.94 


14.47 





4 

46443 

0.72 



14.30'' 

13.86" 




5 

46768 

0.03 

16.91 

15.98 






6 

46801 

0.06 

16.31 

15.69 






6 

46813 

0.08 

16.37 

15.68 






6 

46827 

0.09 

16.28 

15.65 






6 

46835 

0.10 

16.16 

15.52 






6 

46837 

0.10 

16.46 

15.71 






6 

46845 

0.11 

16.19 

15.62 






6 

46873 

0.13 

16.67 

15.81 






6 

46876 

0.14 

16.42 

15.65 






6 

46877 

0.14 

16.19 

15.54 






6 

46895 

0.16 

16.52 

15.70 






6 

46904 

0.16 

16.14 

15.54 






6 

47127'" 

0.38 

15.78-' 


14.54'" 





7 

47556 

0.79 



14.65'" 





7 

49309"' 

0.47 

15.96 

15.31 






8 

49398 

0.54 





12.54 

11.66 

11.33 

9 

50122 

0.25 




13.89 

12.35 



10 

50417 

0.51 




14.00 

12.42 


11.30 

10 

50894 

0.96 





12.51 

11.65 

11.33 

11 

51904 

0.92 





12.57 

11.72 

11.35 

12 

51927 

0.94 





12.54 

11.66 

11.31 

12 

53332 

0.28 





12.41 

11.69 

11.35 
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Notes: mean date, ^ Kron-Cousins magnitude, ^ photographic F magnitude, photographic J magnitude 

References: (1) SuperCOSMOS Sky Surveys (Hambly et al. 2001); (2) Kafatos et al. (1983) (3) USNO-A2.0 Monet et al. (1998); (4) The USNO-Bl.O 
Catalog (Monet, Levine & Canzian 2003); (5) Kontizas et al. (2001); (6) Lawson et al. (1990); (7) The GSC 2.2 Catalogue (STScI 2001); (8) SPM 4.0 
Catalog (Girard et al. 2011); (9) Muerset et al. (1996) ; (10) DENIS Catalogue towai'd Magellanic Clouds (Cioni et al. 2000); (11) 2MASS All-Sky Catalog 
of Point Sources (Cutri et al. 2003); (12) 2MASS 6X Point Source Working Database (Cutri et al. 2006); (13) Kato et al. (2007) 





